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FOREWORD 


This  report  was  prepared  by  Dr.  Robert  H.  Korkegi  and 
Maj  Ronald  A.  Briggs  of  the  Hypersonic  Research  Laboratory,  Aerospace 
Research  Laboratories,  Office  of  Aerospace  Research,  United  States  Air 
Force.  The  investigation  was  carried  out  under  Project  7064  entitled  "High 
Velocity  Fluid  Mechanics.  " 
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ABSTRACT 


Equations  describing  compressible  turbulent  plane  Couette  flow  for 
the  case  of  equal  heat  transfer  at  both  walls  have  been  developed  for  a 
perfect  gas  based  on  the  von  Karman  mixing  length  model.  This  model  was 
selected  because  of  its  good  agreement  with  available  turbulent  plane 
Couette  flow  measurements  of  incompressible  flow  and  furthermore,  lacking 
compressible  Couette  flow  data,  it  is  shown  by  Spalding  and  Chi  that  those 
theories  based  on  the  von  Karman  model  give  best  agreement  with 
compressible  turbulent  boundary  layer  measurements. 

A  sample  calculation  is  given  for  the  heat  rate  to  the  inner  surface  of 
a  rail-guided  slipper  (representative  of  slippers  suppprting  rocket  boosted 
sleds)  based  on  the  assumption  of  a  small  air  gap  between  both  surfaces. 
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NOMENCLATURE 


half-gap  height 

local  skin  friction  coefficient 

local  heat  transfer  coefficient 

constant  empirically  determined 

enthalpy 

eddy  conductivity 

mixing  length 

Eqn.  8 

Mach  number 

Eqn.  8 

pressure 

Prandtl  number 

rate  of  heat  transfer 

gas  constant 

Reynolds  number 

temperature 

velocity 

u/uc 

shear  velocity 
nondimensional  velocity 
distance  from  lower  wall 
nondimensional  distance 


mixing  length  constant 


jj  eddy  viscosity 

T  shear  stress 

Subscripts 
b  half -gap 

c  properties  at  midpoint  of  gap 

d  gap  width 

r  upper  wall 

t  turbulent 

w  lower  wall 

€  eddy 

®  free  stream 


INTRODUCTION 


Turbulent  plane-Couette-like  flow  has  been  observed  experimentally 
1  2 

by  Reichardt  and  Robertson  in  apparatus  utilizing  endless  moving  belts. 

3 

Both  authors,  as  well  as  Squire  ,  present  theoretical  treatments  for  the 
incompressible  case. 

Based  on  experimental  evidence  (cf.  Refs.  2  and  4),  turbulent  plane 

3 

Couette  flow  can  be  expected  to  occur  at  Re<j  =  0  (10  )  where  d  denotes  gap 

2 

height.  In  the  data  of  Couette,  presented  by  Robertson  ,  laminar  flow  was 

3 

observed  for  Re^  =  400  and  turbulent  flow  for  Re^  =  10  between  concentric 
cylinders  with  a  spacing  of  1.  7%  of  the  cylinder  radius. 

Interest  in  compressible  turbulent  plane  Couette  flow  is  stimulated 
by  such  problems  as  the  rail-guided  slippers  supporting  rocket-boosted 
sleds  accelerated  to  speeds  of  several  1000  f.  p.  s-.  at  the  Test  Track 
Directorate  of  AFMDC,  Holloman  AFB.  Fully  developed  shear  flows  can  be 
expected  in  the  air  gap  between  slipper  and  rail  surfaces,  and  the  Reynolds 
numbers  based  on  gap  height  are  large  enough  for  turbulent  flow  to  prevail. 

TURBULENT  FLOW  MODEL 

For  lack  of  experimental  data  on  compressible  turbulent  plane 

Couette  flow,  the  turbulent  law  selected  is  chosen  as  that  which  gives  best 

agreement  with  turbulent  boundary  layer  data  at  supersonic  speeds, 

compatible  with  good  agreement  with  the  low  speed  skin  friction  and  velocity 

2 

profiles  for  turbulent  Couette  flow  measured  by  Robertson  .  The  data  of 
Robertson  were  used  because  of  the  large  number  of  measurements 
presented  and  their  good  agreement  with  the  earlier  work  of  Couette. 

According  to  the  excellent  survey  of  theories  for  turbulent  boundary 
layer  flow  of  Spalding  and  Chi^,  those  based  on  mixing  length  give  best 
agreement  with  experiment  for  both  the  adiabatic  and  nonadiabatic  wall  case. 

7 

Two  of  the  three  theories  giving  best  agreement  are  those  of  Wilson  and 

g 

Van  Driest  ,  both  based  on  the  von  Karman  mixing  length. 


On  the  strength  of  these  findings  coupled  with  the  good  agreement 
with  Robertson1  s  data,  the  von  Karman  mixing  length  is  selected  as 
turbulent  flow  model. 


GOVERNING  EQUATIONS 

For  Couette  flow  the  equations  of  motion  reduce  to 
Momentum  T  =  const  =  Tw  (1) 

Energy  q  -  T  u  =  const  =  qw  (2) 

Where  T  and  q  are  taken  to  denote  the  sum  of  laminar  and  turbulent 
shear  and  heat  transfer  respectively,  the  subscript  w  denotes  values  at  the 
lower  wall,  and  u  denotes  the  time  average  value  of  the  velocity.  A  sketch 
of  the  flow  is  shown  in  Fig.  1. 

Note  that  r  =  const  is  a  simplifying  approximation  in  most  turbulent 
boundary  layer  theories,  whereas  it  is  exact  for  Couette  flow. 

In  the  compressible  analysis,  only  the  symmetrical  case  of  equal 
wall  temperatures  and  heat  rates  will  be  considered,  which  implies  the 
velocity  profile  will  be  antisymmetric  about  the  midpoint  of  flow.  Thus, 
following  Fig.  1  and  utilizing  Eqns.  1  and  2 

^r  ^w  ~  Tw  uc  ^ 

where  uc  =  /4  ur  is  the  midpoint  velocity,  the  subscript  r  denoting  values  at 
the  upper  wji11. 

Neglecting  the  laminar  contribution  to  shear  and  utilizing  the 
von  Karman  mixing  length 


du 


where  K  is  a  mixing  length  constant  whose  value  is  taken  as  0.4  ,  one 
obtains  the  differential  equation 


dfu+ 
dy+ 1 


0 


(4) 


2 


where  y+  and  u*  are  nondimensional  height  and  velocity  defined  as  follows 


(4a) 

uT  =  Vrw/pw  is  the  shear  velocity,  and  Reynolds  number  and  skin 
friction  coefficients  are  based  on  midpoint  values,  b  and  uc  --  see  Fig.  1  -- 
and  wall  viscosity,  and  defined  as  follows 

ucb  Cf  T  w  ^  uT  \  2 

Re,  =  -  ;  —  =  - T  =  ; - - 

b  2  pw  ^c/ 


With  the  boundary  condition  u+  (y+  =  0)  =  0,  the  general  solution  of 
Eqn.  4  yielding  the  velocity  profile  between  the  fixed  wall  and  tlie  midpoint 


of  the  flow  is 


/■»u  *  £=  \  i/ 

y+  =  C  \  exp  (K  \  (£-  ,  /2  d$l  du+ 


If  one  applies  a  second  boundary  condition  of  (du/dy)y  _  q  =  Tw/pw 
at  the  wall,  e.  g.  , 


>y=0 

one  obtains  a  constant  of  integration 


C  =  1. 


However,  as  the  flow  near  the  wall  is  dominated  by  laminar  motion,  one 
cannot  rigorously  apply  a  boundary  condition  based  upon  a  turbulent  shear 
flow  model.  Therefore,  in  accordance  with  standard  turbulent  boundary 
layer  notation  (cf.  Spalding  and  Chi), 


3 


is  used  where  the  value  E  is  determined  empirically.  El  K  may  be  viewed  as 
a  scale  or  stretching  factor  on  the  coordinate  y.  Thus, 


The  density  ratio  p/pw  *or  the  compressible  case  is  related  to  the 
velocity  distribution  by  means  of  the  energy  equation,  Eqn.  2.  Introducing 
an  eddy  viscosity  pe  ,  eddy  conductivity  kg  ,  and  turbulent  Prandtl  number 
pr^  —  cp  fig  /  kg  ,  the  turbulent  shear  and  heat  rate  are  written  respectively 

du  ,  kg  dh  _  Mg  dh 

Tt  ~  ^g  dy  an  qt  ~  cp  dy  P*g  dy 


Neglecting  the  laminar  contribution,  from  Eqn.  1,  Tw  =  pe  (du/dy)  -  const 
hence  q  =  qt  =  -  (Tw/Fr€  )  (dh/du),  and  therefore  Eqn.  2  is  written 


Tw _  dh 

Pr,  du  Tw  U  ~  'qw 


Assuming  Pr£  =1,.  and  introducing  Eqn.  3  for  the  case  of  equal 
heat  transfer,  its  solution  is 


..2 

h  =  hw  +  uc  U,  -  y 


(7? 


A  Mach  number,  M_,  based  on  midpoint  velocity  and  wall  enthalpy  is 


defined 


m*  =  ^  Mc2  = 


(8) 


Equation  7  can  thus  be  written 


where 


— —  =  1  +  m2  (2u  -  u2) 
hw 


_  u 

U  =  - 

U, 


(7a) 


4 


Now,  limiting  the  analyeis  to  the  case  of  a  perfect  gas  -  -  p  = 
pRT  and  noting  that  the  pressure  is  constant  throughout  the  Couette 


£_  = 
Pw 


1  +  m*  {2u 


u  -  u2)  1 


INCOMPRESSIBLE  CASE 


For  p  =  pw  =  const,  the  solution  of  Eqn.  6  is 


and  the  skin  friction  law  (u  =  u_  for  y  =  b) 


-Vi-  =  j? (> +  v-'Jt  r%) 


Best  fit  with  the  experimental  data  of  Robertson  is  obtained  for 

E  =  8.  8  (Note  that  E  varies  from  about  1 1  to  1 3+  in  various  turbulent 

boundary  layer  theories  based  on  mixing  length,  cf.  Ref.  6),  and  further- 

3 

more,  since  turbulent  Couette  flow  has  been  observed  only  for  Re^,  >  10 
for  which  (E  VCf  f  2  Rej,)i  >  200,  Eqn.  11  can  be  approximated  as 


-M‘ 2-51"  \8-8‘Vf  Reb) 


i  2  ,  /Cf  \ 

“\/c7*  5'44+  2-Slnl~VT  Re„) 


This  equation  is  plotted  in  Fig.  2  and  gives  very  good  agreement 
2 

with  Robertson1  s  measured  values  as  well  as  those  of  Couette  given  in  the 
same  reference  over  the  whole  range  of  Re^.  Also  shown  in  Fig.  2  is  the 
laminar  law  Cf  Rej.  =  2.  Note  that  Eqn.  12  is  almost  identical  to  the  skin 


5 


friction  law  given  by  Robertson  based  on  Reichardt1  s  solution  for  the 
velocity  distribution  assuming  a  parabolic  variation  of  eddy  viscosity 

5.5+  5.75  1og10  (ltcb  ) 

Velocity  profiles  from  Eqn.  10  with  K  =  0.  4  and  E  =  8.  8,  are  given 
3  4  5 

in  Fig.  3  for  Reb  =  10  ,  10  ,  and  10  and  compared  with  Robertson' s 

4 

measurements  for  Reb  =  10  . 


DEDUCTIONS  FROM  THE  INCOMPRESSIBLE  LAW 


If  the  velocity  profile  very  close  to  the  wall  is  examined,  the  log 
term  in  Eqn.  10  can  be  expanded  for 


which  has  the  desired  linear  laminar  character  and  precisely  the  laminar 
form  for  E Ik  -  1  except  that  the  velocity  profile  has  a  much  steeper  slope 
as  [  (Cf/2)  Reb]turb  >>  1  compared  with  [  (Cf/  2)  Reb]lam  =  1. 

It  is  thus  found  that  the  von  Karman  mixing  length  model,  which 
depends  only  on  the  local  flow  characteristics,  gives  the  features  of  a 
laminar  sublayer  when  applied  near  a  wall. 

Also,  taking  the  slope  of  the  velocity  profile  from  Eqn.  10 


A  /U  ' 

d  \TT~  I 
\uc  / 

<*&-/ 


6 


At  the  wall 


d 


d 


u  \ 

\*c) 


y  = 


0 


E 

K 


Si 

2 


Re 


-b 


as  from  Eqn. 


13  and  at  the  midpoint  y  = 


d  ( 


/u 


ur 


E 

K 


l1- 

\b  /' 


y  =  b 


/Cf 

1  +  E  “VT  Reb 


It  is  also  interesting  to  note  the  form  of  the  eddy  viscosity  for  the 
incompressible  case  obtained  from  the  slope  expressed  above  and  rewritten 
as  follows 


/  .  E  \  du 

Tw-  P\ji  (‘'w  +  kut  K  y)  IE  \ 

Here  again  appears  the  laminar  form  and  the  constant  E/k  as  a  scale  or 
stretching  factor.  Note  that  the  form  of  the  eddy  kinematic  viscosity  from 
von  Harman1  s  , mixing  length  model  is  linear  for  incompressible  Couette 
flow 


uT=  KuTy 

As  a  matter  of  interest,  the  ratio  of  turbulent  to  laminar  viscosity  is 

turbulent  viscosity  ^UT  y  — ^  /C f;  y 

■; - : - : - =  -  =  E  \J  —  ReK  f- 

laminar  viscosity  Vw  2  D  b 

Thus,  the  approximation  made  in  Eqn.  12  corresponds  to  neglecting 
the  laminar  contribution  to  the  shear. 

The  form  of  the  eddy  viscosity  is  actually  not  too  important  --  what 
is  important  is  the  scale.  Thus,  if  a  parabolic  eddy  viscosity  distribution 
is  taken,,  rather  than  a  linear  variation,  i.  e.  vT  =  KuTy(l-u/2b) 


7 


(following  Reichardt),  the  skin  friction  becomes 


which  differs  from  Eqn.  11  only  by  a  second  order  term. 


COMPRESSIBLE  CASE 

In  a  manner  analogous  with  turbulent  boundary  layer  development, 
it  is  assumed  that  the  incompressible  scale  factor  holds  for  the  compress¬ 
ible  case,  i.  e.  E  =  8.8. 

Introducing  Eqn.  9  for  the  density  ratio  into  Eqn.  6,  one  obtains 


£  v+  -C 


„  exp  {*  I 


U  .  -Vf“+ 

[  1  +  m2  ( 2u 


iu  -  u2)]  ^dul 


duf(14) 


J 


Integration  yields  the  velocity  distribution 
(i  +  ~T  m2 )  f^E  “V  Reb  ^  =  j^'l  +  m2  (2u  -  u2) 

1  /Of  2  ,,  _J  I  fl  1  T  .  - 1  m  .  -1  (1  -u)  m  1 

K  V  2  J  1  vcf  m  L  Nli+rri*  ^1+m2  J 


-1  +  ±  m2 

K  v  2 


(15), 


and  the  skin  friction  law  (u=  1,  y/b  =  1) 

(‘  +  p-  T  m!)  (E  ReW  =  *'TT 


mz  exp 


.'2  1  .  - 1 


/  __  _Z_  1 

(K^  Cf  m 


tan  m 


- 1  +  -  m2 

K  *  2 


(16) 


for  compressible  flow. 


Eqn,  16  can  be  approximated  in  similar  fashion  as  Eqn.  11  and 

written 


m 


1  + 


,  "l 
tan-  m 


In 


1 


Cf 

2 


m 


\/ 1  +  m2 


E 


(16  a) 


Note  that,  as  m  — »  0,  (tan  *m/m)  — *  1,  and  Eqns.  15  and  16  reduce  to  the 
incompressible  form  of  Eqns.  11  and  12  respectively. 

Fig.  4  shows  the  compressible  turbulent  plane  Couette  flow  skin 
friction  law  given  by  Eqn.  16a  as  well  as  the  laminar  one  obtained  assuming 
a  linear  viscosity-temperature  relationship,  and  Pr  =  3/4,  for  equal  wall 
temperature  and  heats  rates 


1  + 


Zzl 

3 


PrMc2 


Note  that,  contrary  to  the  turbulent  case,  the  laminar  skin  friction 
increases  with  Mach  number  due  to  increasing  wall  cooling. 

Also  shown  in  Fig.  4  is  the  heat  transfer  coefficient  whose 
definition  is 


Ch  = 


^w _ 

PWUC(K  -  hb> 


where  h-D  is  the  midpoint  enthalpy.  From  Eqn.  7  --  h^  =  hw  +  uc2/  2  -  -  and 
Eqn.  3,  it  follows  that  C ^  =  Cf.  Alternately,  the  heat  transfer  rate  to  the 
walls  may  be  calculated  directly  from  Eqn.  3. 

9r  —  ~^w  ~  2  ^f^wuc  (17) 


The  ratio  of  compressible  to  incompressible  skin  friction  coefficients 
for  various  Reynolds  numbers  is  shown  in  Fig.  5.  ^ 

Typical  turbulent  velocity  profiles,  given  by  Eqn.  15,  are  shown  in 
Figs.  6  and  7  for  various  Mach  numbers  and  Reynolds  numbers.  The 
density  and  enthalpy  profiles,  obtained  from  Eqns.  9  and  15,  are  shown  in 
Figs.  8  and  9. 
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The  very  small  variation  of  density  across  most  of  the  profiles,  as 
shown  in  Fig.  8,  suggests  that  the  simplifying  approximation  of  assuming 
constant  density  at  the  midpoint  value,  p^,  in  the  integration  of  Eqn.  6  might 
be  useful. 

Thus,  with(p^/pw)=  (1  +  m2)  *  from  Eqn.  9,  Eqn.  6  yields 
essentially  the  incompressible  velocity  profile  modified  for  compressibility 


—  =  i_V?"  In 

uc  K  v  2 

and  an  approximate  skin  friction  law 


E~VT  Reb 

\/ 1  +  m2 


(18a) 


(18b) 


The  velocity  profiles  given  by  Eqn.  18a  are  somewhat  fuller  than 
those  given  by  Eqn.  15,  with  increasing  departure  at  higher  Mach  numbers. 
As  for  the  skin  friction  law,  C£,  given  by  Eqn.  18b,  is  lower  than  that  of 
Eqn.  16a  by  about  5%  at  Mc  =  1,  15%  at  Mc  =  2,  and  increases  to  about  40% 
at  Mc  =  8  for  103  g  Refa  g  106. 


Example : 

Consider  a  slipper  bearing  riding  along  a  rail  at  supersonic  speed 
with  a  small  air  gap  between  the  slipper  and  rail  surfaces  so  that,  in 
cross  section,  the  flow  may  be  considered  two  dimensional  in  the  first 
approximation.  Assume  that  at  some  distance  behind  the  leading  edge  of 
the  slipper,  the  flow  in  the  gap  becomes  fully  viscous  and  Couette  flow 
develops.  A  sketch  of  the  longitudinal  cross  section  is  shown  in  Fig.  10 
indicating  the  approximate  bow  shock  of  the  slipper. 

For  purposes  of  illustration  it  is  assumed  that  the  air  entering  the. 
gap  traverses  a  normal  shock  portion  of  the  bow  shock,  and  the  pressure 
in  the  gap  is  constant  and  corresponds  to  that  behind  a  normal  shock. 
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Assume  u  =  4500  £.  p.  s.  so  that  u  =  u-  /  2  =  2250  f.  p.  s.  ,  and 
the  gap  height  is  1  /  8"  so  that  b  =  1/16".  Taking  atmospheric  temperature 

as  70  *F,  Tffl  =  530  #R;  the  speed  of  sound,  a=  N/yR  =  1128  f.  p.  s.  ,  so 
=  4.  0  and  =  2. 

From  the  normal  shock  relation  for  y  =  1.4 

Pgas  .  -1 

P  6 

*co 

Here,  assuming  atmospheric  wall  temperature  for  slipper  and  rail  -- 
Tw  =  =  530°R  --  the  wall  density  is 


.  tlSS.  „  A-  (28Mc  ;1  \  -  0; 


w  R  T 


slup 

0430 


The  Reynolds  number  based  on  b  is 

Pw  uc  b  / 

Re,  =  —  =  1.  38x10° 

b  Mw 

so  that  the  Couette  flow  can  be  expected  to  be  turbulent. 

For  Mc  =  2  and  Reb  =  1.  38x1 0b,  the  turbulent  skin  friction 
coefficient,  as  given  by  Fig.  4,  is 

Cf  =  .0014 

The  wall  heat  transfer  rate,  from  Eqn.  17,  is  then 

=  7  r.rPwuc3 


=  450 


ft  sec 


This  heat  rate  is  comparable  to  those  experienced  at  the  stagnation  point  of 
some  bodies  during  reentry.  For  example,  the  stagnation  heat  rates  to  a 
Mercury  capsuie  are  approximately  700  Btu/ ft2  sec. 
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Theoretical  Skin  Friction  and  Heat  Transfer  Coefficients 
for  Turbulent  and  Laminar  Plane  Couette  Flow  Compared 
with  Experimental  Data  of  Ref.  2 


Theoretical  Skin  Friction  and  Heat  Transfer  Coefficients 
for  Turbulent  Compressible  Plane  Couette  Flow  with 
Equal  Wall  Heating 


Compressible  to  Incompres 
nts  for  Various  Re^ 


sible  Skin  Friction 


Fig.  10.  Sketch  of  Slipper  in  Proximity  to  Rail 
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ABSTRACT 


'-Equations  describing  compressible  turbulent  plane  Couette  flow  for  the  case 
of  equal  heat  transfer  at  both  walls  have  been  developed  for  a  perfect  gas  based  on 
the  von  Karman  mixing  length  model.  This  model  was  selected  because  of  its  good 
agreement  with  available  turbulent  plane  Couette  flow  measurements  of  incom¬ 
pressible  flow  and  furthermore,  lacking  compressible  Couette  flow  data,  it  is  shown 
by  Spalding  and  Chi  that  those  theories  based  on  the  von  Karman  model  give  best 
agreement  with  compressible  turbulent  boundary  layer  measurements. 

A  sample  calculation  is  given  for  the  heat  rate  to  the  inner  surface  of  a 
rail-guided  slipper  (representative  of  slippers  supporting  rocket  boosted  sleds) 
based  on  the  assumption  of  a  small  air  gap  between  both  surfaces.  '  > 
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